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Q^' Abstract. In the framework of Wilson Chiral Perturbation Theory T], we study the 

^ ■ effect induced by a twisted Wilson term, as it appears in Twisted Mass QCD ,2| (with 2 

degenerate quarks). In particular we consider the vacuum orientation and the pion masses. 
The computations are done to NLO both in the mass and in the lattice spacing (i.e. to 
^ \ 0{a?)). There are no restrictions on the relative size of lattice artifacts with respect to the 

' physical mass, thus allowing, in principle, to bridge between the physical regime and the 

unphysical one, where lattice artifacts tend to dominate. The inclusion of 0{a?) lattice 
artifacts can account for the splitting of degeneracy of the three pion masses. Moreover 
O(a^) terms are necessary to model non trivial behaviors of the vacuum orientation such as 
possible Aoki phases. It turns out that these last two phenomena are determined by the 

same constant. 
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1 Introduction and Motivations 



Twisted mass QCD (tmQCD) [2j is a promising regularization for Lattice QCD, because it 
might succeed in improving on many of the most annoying problems of Lattice QCD with 
Wilson-type fermions. It should be able to provide an easy way to reduce lattice artifacts, 
allow simulations of lighter quark masses and even of non degenerate quark masses, still with 
a positive determinant [31 E] • This is accomplished, basically, by introducing a new degree 
of freedom in Wilson Lattice QCD: the chiral-flavor orientation Q of the subtracted Wilson 
term^: Wcr = — ^V* ■ V + Mcr{g)- The most general Lattice Action, that we refer to, is: 

S^[U, ^, g; Mg] = Sg[U; g'] + ^ V'x [7 ■ V + e^^^'^'^X. + M,],, (1) 

xy 

Wilson Chiral Perturbation Theory (WChPT) PI E] has proved to be a crucial tool for 
the analysis of lattice data [El 13 El E] • WChPT has already been extended to the tmQCD 
case in [IHIIIII, including 0(a) terms. 

There are good reasons to consider also O(a^) corrections. In fact, as proposed in j^l 
E], the special choice of twisting u = 7t/2 leads to automatic 0(a) improvement of most 
observables, and thus the leading lattice artifacts are necessarily of O(a^). Even if a; 7^ 7r/2, 
the three pions of {Nf = 2) tmQCD are still degenerate to 0(a). However tmQCD breaks 
explicitly the flavor symmetry, for which the pion mass splitting is an important signal. 
Such splitting appears only at O(a^). Finally the inclusion of O(a^) terms in the Chiral 
Lagrangian is necessary in order to describe the phase diagram of Wilson Lattice QCD at 
very small masses ^21 Ej. It is natural to ask whether this is crucial also in presence of 
twisting, and what is the scenario in this case. 

We have a problem, however, with Twisted Wilson ChPT. In order to fully profit of an 
expansion to O(a^), we would like to treat the mass m and the lattice spacing a on the 
same level, without requiring m » aA^ (A will be defined later). But if we do that, we 
do not have anymore a fixed vacuum around which we can expand^. In fact when m and 
a point to different directions, the limit m — >■ and a — do not commute. In principle, 
of course, the continuum limit should be performed first, and the chiral limit afterwards. 
But the whole idea of WChPT is to reproduce the practical situation, that occurs in Monte 
Carlo simulations, where both quantities are never small enough. Therefore we want to keep, 
as much as possible, a general set up. As we will see in section El this is just a technical 
difficulty that can be easily overcome. 

The purpose of this work is to provide a support to lattice unquenched simulations of 
tmQCD which have been already started E]. For a recent review of ChPT at finite a 
in general, we refer to JH]- Other works have been presented very recently studying related 
problems f71^1Ell2ni- We hope that from so much activity a clearer picture can emerge. 

^We use a standard notation: a is the lattice spacing, r the Wilson parameter, g the gauge coupling, 
Mcrig) the critical mass, V (V*, V) is the forward (backward, symmetrized) covariant derivative, the 
Pauli matrices. 

^We adopt the convention of calling "vacuum" the saddle point of the Chiral Lagrangian around which 
the perturbative expansion is performed. 
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In section 121 we write explicitly the Chiral Lagrangian that we use, and give some com- 
ments. In section 121 we compute the vacuum orientation, and study the special behavior near 
the critical region. In section H] we compute the pion masses and in section El we give a final 
discussion. 

2 Chiral Effective Lagrangian 

The low energy description for Lattice tmQCD that we use in the present work, is given by 
the effective Lagrangian C^ = Clo + ^^nlo lUl HI 1^1^: 

Clo = ^ (S^SS^St) + ^ (mSt + Sm) + ^ (aS^ + Sat) , (2) 



+ L7 (mE^ - Em)^ + W7 (mT) - Em) (aT) - Ea^) + laT) - Ea^)^ + 



The field E(x) is related to the pion field by 

i T-7v(x) 

E(x) = e ~, 

where F ~ 93MeV is the pion decay constant, and tj are the Pauli matrices. The coefficients 
Li [2Ij are the usual dimensionless Gasser-Leutwyler's Low Energy Constants (LEG). The 
other dimensionless LEC's Wi and the Wl - introduced in [H E] - describe the effect of 
the lattice artifacts which appear in the Wilson Lattice formulation of QCD (WLQCD) 
[22] • Space time indices are summed according to the Minkowski metric (the Euclidean 
Lagrangian can be recovered by introducing a minus sign for each (9E(9E) term). In order 
to describe the tmQCD case, the parameters rh and a are 2 by 2 matrices in ffavor space. 
In this paper we always adopt the "physical basis" of tmQCD (see |2lj), in which case^ 



/m \ . / a{cos uj + tsm uj) \ , . 

m = 2Bo , a = 2Wo ^ ' , .... (3) 

\ m / \ a(cosc<; — z smcj) / 

For later convenience we introduce the dimensionless parameters, (as already done in [Zj) 

_ 25om _ 2Woa _ p 



^We also follow the choice of where explicit reference to csw has been dropped. 
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We will never use the "twisted basis", in which lattice artifacts are oc I2. However the 
conversion is straightforward: /i, = msin(a;) and an ordinary mass nig = mcos{uj). 

The Lagrangian above represents an expansion for small momenta p^, small quark masses 
m and small lattice spacing a. In particular we expect that it gives a good description of 
WLQCD when p^, Bom, Woa << A^, where A ~ 1 GeV. In Q all the terms necessary 
for a calculation to NLO of pion masses are included. We refer to ^1 El for a full 
justification why no other terms need to be included at this order. As usual we choose 
-Bo > 0. In principle we do not know the sign of Wq, but it can always be made positive by 
a redefinition of — + vr. Since we are going to study the whole range in u, the choice 
Wq > is not a real limitation. For the same reason, we can restrict to a > 0. Finally the 
system at negative m has a mirror description in the system with positive m. Thus in the 
following we will assume a, m, t] > 0. 

Here m corresponds to "some" definition of the renormalized quark mass, and a to "some" 
definition of the lattice spacing. By construction of the Chiral Effective Lagrangian the 
leg's Li, Wi, Wl do not depend on the size of m or a, but they do depend on the particular 
definition which is chosen for m and a (for instance based on pion mass rather than PCAC)^. 
In the continuum limit the physical Lj are expected to loose also this dependence, but we 
have no reason to think that the same will happen for the W^s. For instance if we redefine 
the mass and lattice spacing: 

m* = m + d, d* = d, (4) 
we find that the LEC's are correspondingly changed in 

-^4,5,6,7,8 = -^4,5,6,7,8, (5) 

= 1^4,5-^^4,5, 
^6*7,8 ~ W^6,7,8 — 2L6j^8, 
^'6,7,8 = W^'6,7,8 — 116,7,8 + -^6,7,8, 

and the LO 0(a) terms disappear, leaving only terms of order 0{p^ a) and 0{ma). We 
will occasionally employ m* later, which has both technical advantages (when computing 
logarithms) and provides an interesting point of view. But we mainly use the parameters 
m and a, as in 0, because they are in the physical basis, and have otherwise a generic 
form that represents your preferred definition of m and a (on the lattice m* is not always 
accessible). 

To avoid confusion, let me state clearly that an 0(a) redefinition of the mass as in 
does not reshuffle the Symanzik expansion of course, which is well defined (given a lattice 
action), since the mass term and the Pauli term are different operators. It is only at the level 
of LO Chiral Lagrangian that the ambiguity appears, since at that level we have essentially 
only one independent operator. 

In order to make the contact with lattice simulations as clear as possible, let me add 
a last remark (see pQ for a slightly different formulation of the same concept). The low 

''Of course the Wi and VF/ depend also on the choice of the lattice action. 
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energy representation Q of Lattice QCD is meant to hold for the renormahzed Lagrangian. 
Changing the lattice spacing a (typically through a change of /3, or the Wilson parameter r) 
has big effects, at the level of the cutoff scale, which cannot certainly be taken into account 
by the expansion in However this possibly strong dependence is compensated by the 
dependence on a in the renormalization constants, delivering a renormalized action whose 
"residual" a dependence is parameterized in 

3 Vacuum orientation 

The first problem, that we have to consider, is the determination of the minimum of the 
potential derived from when the field is set to a constant 2 by 2 unitary matrix Eq 
and the derivatives are set to zero. In a typical ChPT calculation the vacuum orientation 
is changed by the introduction of a source. However such sources are always extrapolated 
to zero in the final result, because they are needed only to derive observables. Even non 
degenerate quark masses do not produce a rotation of the vacuum in the chiral limit, since 
vector-flavor symmetry is not spontaneously broken. As a result Sq = 1 in the final expres- 
sions. The case of ChPT for tmQCD is peculiar, because a vacuum rotation is produced by 
a non zero a and uo, which are not sources, but true parameters of the model that we want to 
describe. Moreover the vacuum orientation Sq does not go to 1 even in the limit a, m ^ 0. 
In fact Eq depends on the ratio a/m m which the limit is attained. This is obvious from 
the physical point of view, because the direction in which chiral symmetry is spontaneously 
broken, depends on the direction, in chiral-fiavor space, where the mass term is pointing. At 
leading order in the effective theory the twisted lattice artifacts are nothing but an effective 
mass term. In other words the limit a — and m — do not commute, and a double 
expansion in a and m cannot be defined around a single vacuum. 

Of course in the physical regime the continuum limit a — should always be taken 
before the chiral limit m — > 0. Here the assumption r] « 1 makes sense, and the point 
of minimum Sg deviates from 1 only by small corrections of 0(a). This is the approach 
adopted in piU lllj. However, in this paper, we want to study also the regime where lattice 
artifacts are not necessarily much smaller than the mass term. This is because we want to 
represent the practical situation that occurs in Lattice Monte Carlo simulations, where the 
continuum limit and the chiral limit cannot be completely and safely separated. In fact, this 
is the whole idea of WChPT. Moreover we also would like to describe the deep unphysical 
regime where lattice artifacts tend to dominate over the mass term, which can have very 
interesting features, as first noticed in [T^ . 

The strategy of the present computation is to fix a ratio r] first, compute the corresponding 
vacuum So(?7), and perform an expansion in x for each fixed rj. This is nothing but a 
convenient organization of the computation, which let us read off, at the end, the results for 
different regimes. Equivalently one could perform a (?7-dependent) change of variables to m* 
in (jU and expand around the fixed vacuum determined by m* (which is however not fixed 
in the physical basis). 
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Figure 1: LO solution 6lo as function of a;. 



In order to find the minimum of the potential we use the parameterization: 

Eq = I2 cos(6') + ir^ sin(6') cos(0) + iri sm{9) sin(0), (6) 

Having introduced a twisted term in the direction, we expect a vacuum rotation in the 
same direction. However we cannot - a priori - exclude a more complicated pattern. Since 
the directions ri^2 are fully equivalent, the parameterization (jHl) is completely general. 



3.1 Potential at LO 

At LO the potential is (proportional to): 

'^xLO = —V cos(u;) cos(6') — cos(6') — t] cos(0) sin(cc;) sin(6') (7) 
which has a minimum for 6 = and 



T] cos(ti;) + 1 



if 1 + ?7cos(cj) > 

if 1 +i]cos{uj) < && sini 
if 1 +r]cos{iu) < && sin(u;) < 



Olo = tan"^ ( ^ ^^^Y 1 ) + ( ^ if ^ +r]cos{uj) < && sin(cj) > 



-TT 



In our convention the Image of tan~^ = (— 7r/2, 7r/2). The trivial orientation {6 = 0) is 
recovered whenever ci; = 0,2'7ror?7<<l. When 77 is not negligible, the vacuum undergoes 
a fluctuation as in figure [T] (left), but it is never rotated more than 7i/2, as long as 77 < 1. 
When the lattice artifacts dominate rj > 1, the vacuum undergoes a full rotation, tuned by 
uj, as in figure ^ (right). We see that for t] = 1 and u) = n the system is critical: it is the 
point where the leading order "effective mass" m* changes sign. Here the LO terms in the 
potential cancel, and the NLO terms become relevant. Around this point one expects that 
phenomena like ^21 EI niay show up. The fact that this happens at = vr instead of = 0, 
is simply a convention related to our choice Wq > 0. 



6 



The stationary point (jH)) is indeed a minimum, since the second derivatives are always 
positive: 

—V^ =^il^ + 2cos{uj)r^ + l (9) 



h2 



^772+2 cos(a;)r7+l 



Notice that ^xi^) ~ ^ each. 6 (this will be true also at NLO). Therefore the mixed 

second derivatives are always zero in = 0, tt, which simplifies considerably the study of the 
stability. Since only one maximum and one minimum are possible in the function ((7j), this 
conclude the study of the LO vacuum. 



3.2 Potential at NLO 

If we add the NLO terms to the potential, we have: 

= V^LO - x[4 cos(2 9) LsG + 4 (cos(cj) cos(2 6) + cos(0) sin(cj) sin(2 6)) Wse t] 
+ (cos(2^) + cos(2c<;) + 3 cos(2c<;) cos(2^) + 
+4 cos(2(/))sin2(tu)sin2(e) + 4 cos(0) sin(2cj) sin(2^))iy^6^2] ^^o) 

where we have defined XgQ = Xg + 2X6, for X = L,W, W, which are the only relevant 
LEG'S, here. 

First of all we deal with the possibility of a non trivial stationary point in </>. One finds 
that, besides = 0, vr, other possible solutions of = are: 

_i / 1 + 8x cosi9)Ws6 + 16VX cos(cu) cos(^)iy'86 



However, in the regime where ChPT is applicable, the relation xX << 1 and XV-^ << 1 
must hold for any LEG X. Therefore (jllj) can never be a solution for our problem. It is 
known ^21 ^| (and we will see it later in the context of tmQGD), that NLO terms can 
actually produce new solutions for the minimum of the potential, but only if LO terms are 
subject to cancellations, which are not possible in (fTTj). This seems to exclude the possibility 
of a phase of spontaneously broken fiavor symmetry (i.e. along a direction different from 
the one of twisting, when twisting is non zero), as it has been very recently suggested |17J. 
Although the general functional form of (llUj). would allow it, it is not within the reach of 
our GhPT representation. This is an interesting question, which certainly deserve further 
numerical studies, and may find a description within an alternative analytical approach. 
However, in this paper, we set in the following = (the choice = tt is equivalent, if one 
correspondingly changes 6 —0). 

The NLO correction to the solution (jSI) is: 

= Olo + ^Xiv' + 2 cos(cu)r^ + 1)'^ [ - sin (2tan-^( ^ ""'""^^^ ))L,, + (12) 

7] cos(ct;) + 1 



7 



+ sin (2a; - 2tan-^( ''^' ''.' Wsev' + sin (a; - 2tan-^( ''^'P'^' Wsev] 

In figure|21we plot (in dashed-blue) the solutions (fT^ for some typical value of the parameters. 
The parameter 77 is chosen in order to put in evidence the effect of lattice artifacts (77 = 
1/3,3/2), while x ^^nd the LEC's are taken in the range that one can expect from typical 
lattice simulations: (x-^86; XV^se, Xv'^^m ~ {O^^^l/^O}; see the discussion at the end of 
this section). For comparison, next to the solutions ()12|). we plot in full-black the LO solution 
(jHl), and in dotted- red the solution obtained by numerically minimizing the potential (jiup . 
The numerical minimization is performed with the Mathematica routine FindMinimum. In 
particular the comparison between the dotted-red and the dashed-blue curves provides an 
interesting visualization of what one could expect from NNLO corrections. Another insight 
is given by figure El in which the same NLO vacuum solutions 6 are plotted as a function of t] 
for two fixed u = tt/2, 3tt/A. We notice that at = 7c/2 the two NLO descriptions are quite 
consistent over the whole range of rj. The stability of the solution (jl2j) is still guaranteed 
by the positivity of Q. The condition of stability in ^ - including the NLO contributions is: 

— =^v' + 2 cos(a;)r/ + 1 + ^.^,ZUn+i + 2 cos(a;)r^ + l)Lse (13) 

+r7((cos(u;)?7^ + 2t] + cos(u;))py86 + viv^ + 2 cos{uj)r] + cos(2co'))iygg)]. 

However the LO term go to zero when rj 1 and uj —>■ n at the same time. This is the point 
where something new can happen, depending on NLO terms. This is the subject of the next 
paragraph. 

3.3 Critical region 

As already justified we can limit ourselves to consider the potential at = 0: 

V^ = - cos{e) - cos(tu - e)ri - 4x[cos(2e)L86 + cos(cj - 2e)W86V + cos(2(cj - 9))W^Qri'^] (14) 

Notice that NLO terms have a periodicity in 6 which is double as fast as the one in the 
LO terms (i.e. they have period tt instead of 27r). This means that the NLO terms can 
introduce at most a new pair of stationary points (in this case a new minimum and a new 
maximum). In more formal terms, the condition ^V^lo = is a polynomial of degree 2 in 
the complex variable z = e*^, while the same equation for (jl4p has degree 4. The appearance 
of new stationary points is actually only possible if the NLO terms are strong enough to 
compensate the LO term in the stability condition (fT^ . Since we always assume that 
{X = any LEG) are small, the only possibility is when (r^^ + 2?7 cos(c(j) + 1) is small, i.e. when 
both r] 1 and u —>■ n. 

In order to describe this region we choose a convenient parameterization: 

77 = 1 + xSt], w = tt + xSuj, (15) 
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1 , 1 

e [ll=J, L8i=0, 1 1*86 = 0, I) W'!6 = -2Q-1 



, 1 . 1 2 

e [ll=-J, 186=0, 1] >*!6="-2q"' 1 "'e 



,1 1 .2 

e [ll=J, l8S=2j-: 1 1*86=0, rj W's 




1 , 1 

e [l=J, 186=0, I) 1*86=0, I| W'86=--2j-l 



, 1 . 1 2 

e [i]=j, 18S=0, 1 1*8S = -2Q-' 1 l*'s 



,1 1 .2 





Figure 2: NLO solution for the vacuum orientation 6' as a function of u for various choices 
of the parameters (dashed-blue). NLO solution obtained by numerical minimization of the 
NLO potential (dotted-red). LO solution (full-black). 
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Figure 3: NLO solution for the vacuum orientation ^ as a function of rj for various choices 
of the parameters. Dashed-blue, dotted-red, full-black as before. 
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and we assume that Sui and St] are dimensionless number of the same order of magnitude as 
the LEC's. Roughly speaking one can think of Sr] as a very small ordinary mass, and 6uj as 
a very small twisted mass. In these variables the potential (at first order in x) becomes: 

Vjx = Sr]cos{e) - 8(L86 - Wse + W^,) cos{ey + 6ujsm{e) + 0{x) 

If 6uj = 0, we recover the familiar potential first studied in [13j. This identifies the coefficients 
Ci, C2 and e introduced in [12] with 



C2 = -8{Lse-Wse + W'se), (16) 
-ci = 6r] = F^ 



^2 Woa — Bom 



2(5om)2 
ci/(2c2) = F 



2 Woa — Bom 



32{Bomy[{Lse-Wse + W'se)] 

Of course an overall (positive) factor in ci and C2 is irrelevant. 

It is interesting to see that, although some LEC's depend on the definition of the mass, 
the picture above is stable under such redefinition. We mentioned in section |2l that one can 
redefine the mass as in (j3]) and go to very small m* ~ a^. In this setup one expects that only 
the 0(a2) LEC's W contribute (see [IB]). However from © one finds (Lge - Wsq + W'se) = 
W'Iq, showing perfect consistency between the two descriptions. We can also observe that 
the combination (Lgg — Wse + W'se) is invariant under renormalization (see section The 
constant W'Iq, instead, does not renormalize because there is no LO term in a* which 
produces divergences which need to be subtracted. 

Before introducing a non zero 6uj, we recall the results of the analysis of [TB] : 

1. if Ci > (positive mass m*) and 2c2 < Ci (i.e. [Lge — WgQ + Wge] positive or small 
negative), the solution is 6* = 0. 

2. if ci > (positive mass m*) and Ci < 2c2 (i.e. [Lgg — Wse + W'se] large negative), the 
solution is cos(^) = e > (the chiral condensate has a positive component like in 1.). 

3. if ci < (negative mass m*) and 2c2 < — ci (i.e. [Lge — Wse + W'se] positive or small 
negative), the solution is 6^ = vr. 

4. if ci < (negative mass m*) and — ci < 2c2 (i.e. [Lge — Wse + W'sq] large negative), 
the solution is cos(^) = e < (the chiral condensate has a negative component like in 
3.). 

If we switch on 6uj, the potential does not depend anymore only on cos(6'), but also on 
sin(6'). This breaks explicitly the symmetry ±6. It is clear that the sign of 6 will be opposite 
of the sign of 6uj in order to produce in both cases a negative contribution to the potential 
6ujsm{6). We will not try to find the full solution in this case. In ^| the first correction has 
been computed. Here instead we repeat in figure El a plot analogous to the one in figure |2l but 
now for 1] very near 1. The analytical NLO solution (dashed-blue) display clear instabilities. 
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as expected, being simply a correction to the LO vacuum. It is not anymore relyable, when 
a more comphcated structure of minima can develop. The numerical solution (dotted-red) 
is very interesting: it shows either a jump at a; = vr (see the plots in the second and fourth 
row in figure |2 here C2 > 0), or a more gradual change (first and third row in figured here 
C2 < 0) that involves a temporary transition through a stable minima in ^ = 0. 

The first pattern (a sudden jump) corresponds to C2 > 0. In fact for these parameters 
an Aoki phase is possible, therefore 6' = is a local maximum, and the vacuum can only 
choose between the 2 Aoki vacua at ^ = ±cos~^e. Which one is chosen, depends on the 
sign of o; — TT ~ 5a;. In figure El we follow what happens in details by looking at the potential 
as a function of 9. Going from thick-blue to thin-red, uj goes from tt/2 to vr. The minimum 
generated by the twisted term at = 7r/2 evolves smoothly into the right Aoki minimum. 
At u; = vr the potential is the one typical of the Aoki phase, with 2 minima corresponding to 
9 = ±cos~^(e). If 6uj changes sign, the vacuum tunnels between the two minima, and flips 
the sign of the condensate (ip'y^Tsip) . 

In the second scenario C2 < and no Aoki phase is possible, thus 6' = is a local 
minimum at = vr. When the twisted term is still strong, for instance around u = 7r/2, 
it imposes a vacuum far from 9 = 0. But as soon as the twisted term is weak enough, the 
normal vacuum ^ = becomes again the true solution. In this scenario the transition can 
be smooth, but it can also happen in two separate jumps, signaling a possible coexistence 
of two local minima for some finite value of vr. This is a new phenomena, which does not 
occur in absence of twisting, and it is worth analyzing more closely. The first possibility is 
displayed in figure IHl the minimum at a; = 7r/2 is essentially the LO minimum dictated by 
the twisted term. Approaching u; = vr such minimum is smoothly deformed to the trivial 
one ^ = 0. Although another minimum develops for larger u, it plays no role. The second 
picture can be seen in figure [7| Here the unique minimum at a; = 7r/2 is smoothly deformed 
to the non trivial minimum near 9 = it, which is however unstable, and the transition to the 
true one at 9 = occurs with a jump. In this scenario decreasing the twisted term may have 
the counterintuitive effect of producing a larger (ip'j^Ts'ip) , although only for a short range. 
Unfortunately we could not find a full characterization in terms of LEC's, of these last two 
behaviors. 

Now we should add some comment about the possible values of the parameters that may 
occur in practice. According - for instance - to simulations performed in [71 |H] pion masses 
between 375 and 670 MeV correspond to x between 11 and 36. In the same reference also 
the combination 2Lq + Lg have been estimated^ to be approximately ~ 4.2(1.1) ■ 10~^ (at 
the renormalization scale n = 6F). This combination depends on the renormalization scale, 
and at present no estimate of any W's is available. However this justify our choice (for 
illustration purposes) of plotting the potential in the range x^se = = L, W, W). 

^From data at f3 — 5.1 and no continuum extrapolation has been attempted. 
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Figure 4: NLO solution for the vacuum orientation ^ as a function of uj for various choices 
of the parameters and r/ ~ 1. Dashed-blue, dotted- red, full-black as before. Instabilities are 
clearly displayed in the NLO analytical solution. 



13 



Pot (6) 




Figure 5: Potential V^(6') evolving from lo = -njl (thick-blue) to = vr (thin-red). The 
other parameters are r] = 9/10, xv'^^m ~ ~l/20, L^q = Wse = 0. 



Pot(e) 




Figure 6: Potential Vxi^) evolving from u = 7r/2 (thick-blue) to u = ir (thin-red). The other 
parameters are r] = 9/10, = 1/20, W^q = Wse = 0. 
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Pot (6) 




Figure 7: Potential Vxi^) evolving from uj = 7r/2 (thick-blue) to a; = tt (thin- red). The other 
parameters are 1] = 9/10, xv^^m = 1/20, Lgg = Wge = 0. 

4 Pion masses and pion mass splitting 

In this section we compute the pion masses from the Lagrangian 0, including O(a^) lattice 
artifacts in a twisted direction. The expansion is performed around the vacuum computed 
in section 13.11 and 13.21 As discussed in the previous section, we expect this description to 
be good for any t] and u, as long as x ^i-nd XV ^ire small enough (i.e. we are in the ChPT 
regime), and the combination (ry^ + 2?7cos(co') + 1) is not too small. The latter condition 
means that this expansion cannot describe the "critical" region where both Su and 61] are 
~ C2. There it is necessary to make sure first that the true vacuum has been recognized (see 
section 13. 3|) . 

With these restrictions in mind, we proceed to the computation of the pion mass. For 
that we need only the Lagrangian and we can follow the simple method described for instance 
in [23 (sec. VI — 2). We expand the LO Lagrangian Clq around the vacuum computed 
to NLO order, while for the NLO part of the Lagrangian Cnlo we only need the vacuum 
at LO. The NLO Vacuum - inserted in C^o ^ is necessary to cancel the terms linear in the 
fields 7r(x), coming from the NLO Lagrangian Cnlo- To these two components we must 
add the effective Lagrangian coming from the loop integration Cioop- The computation of 
Cioop is easy, if one observes that the LO Lagrangian expanded around the vacuum at LO 
is equivalent to the untwisted Chiral Lagrangian at LO, without Wilson term, but with a 
mass m* such that: 



This is the new version of the "rule" (x ^ X + p) suggested in |lj to compute the logarithmic 
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term in WChPT from ordinary ChPT. We have then: 

loop 2 ^' ^ 247r2F2 ^ 2 ^ 2A'k^F^ ^ " 

where, in dimensional regularization, we defined, 

A(2Som*) = log(^^^) + ^ - log(47r) +7^-1- 

The pion wave function renormahzation - which is defined by 7r/j(x) = Z~^I'^'k{x) - results: 

S((?7 cos(cj) + 1)^45 + riiji + cos(cc')) W45) 



77^ + 2 cos(ci;)?7 + 1 



^^^^ A(25om*) + O(x^), 

which converges to the result in in the limit // << 1. No new LEG is involved here. If 
we call the pion associated to the untwisted directions t\^2 we find: 

2 



^^^2 + 2cosHr7 + l-x'[^— ^ T^— T (17) 

^ r]"^ + 2 cos(ti;)r7 + 1 

(8(-2L86(r/cos(cj) + 1)^ + (cos(2cj)r/2 + 2r/2 + {jf + 3) cos(cj)r7 + \)L^^ + 

+r/(r7 + cosH)((77^ + 2cos(cu)r7 + 1)1^45 - 2(r/cos(t<;) + 1)1^86 - 2r7(r/ + co^{u)))W'^^)))\ + 

^,(77^ + 2 cosfcu)?] + 1), , / rr; — \ ^, o, 

+^ 32^^^^ ~^ °^ V + + ) + 

Once again in the limit — the result in JUI is recovered. Here we introduced the 
renormalized LEG at the scale as defined in jB]: 



-^45 


— i^45 ~ 


(C5 + 2C4)A(F2) 




= W^45- 


p5 + 2D4)A(F2 


-^86 


= -^86 — 


(Cs + 2C6)A(F2) 


W^86 


= W^86- 


p8 + 2D6)A(F2 


W^86 


= W^86- 


-(C8 + 2C6)A(F2) 


C4 


1 

2567r2 


r - ^ 



Notice that the LEG W[ are renormalized with the same coefficients Cj as the Lj. Notice 
also that, although the functional forms (in r] and a;) in front of the various LEG's are very 
different, the coefficients Ci and Di are such that they sum up to precisely the functional 
form in front of the loop term. This is a non trivial check of consistency for the calculation. 

The most interesting aspect of the formula for the pion mass ()17|) is that each of the LEG's 
enters the expression with a different functional dependence on ri and uj. These are simulation 
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parameters that can be in principle freely changed in Monte Carlo simulations. This allows in 
principle a separate determination of each of the LEG. This somehow reminds the successful 
idea of selecting special combinations of LEG from Partially Quenched simulations [2S1 • But 
the situation is very different here, and whether this could be really exploited in numerical 
simulations, cannot be answered here. 

A twisted Wilson term breaks explicitly the flavor symmetry. Thus one expects a splitting 
in the degeneracy of the pion masses. If we call vr" the pion associated to the twisted 
directions we find: 

ml± - ml, ^ 16r7^x'sin^(^)(^86 - W^86 + W^&) 

F2 + 2 cos(cu)r7 + 1 ^ ' 

^ 64a sm(co') — — — 

The most interesting aspect of this formula is that the pion mass splitting is completely 
determined (to this order) by precisely the same combination of LEG's that decides about 
the possibility of an Aoki phase. In particular m^± > m^o predicts C2 < and no Aoki 
phase. Recent unquenched simulations [131 CS] suggest that there is indeed no Aoki phase 
at any interesting j3. It would be nice to confirm this picture through the measurement of 
the pion masses. Unfortunately m^o is a bit more difficult to measure than the other two in 
tmQGD, since it may have stronger contaminations from heavier states, but it should not 
be impossible. 

As already mentioned, the determinations of 2Lq + Lg, which are available both from 
the lattice |H| and from experiments [2ZI, are not enough for our purposes, because such 
combination alone depends on the renormalization scale, while the combination in (jl8p does 
not. 



4.1 Pion masses near the critical region 

Our final task is to try to make connection between the pion masses computed in fll7ll8|) 
and the computations in (TH]. As already said, we do not expect this to succeed, because 
the vacuum around which we expand is blind to the appearance of a complicated structure 
of minima. However some partial result can be obtained. 

We employ again the expansion (jl5|) that we used before for the critical region. In these 
variables the "charged" pion masses become 

-j^=X [ + +5.y (19) 

while the "neutral" pion is: 



If we recall the identification ()16p and set 5c<j ^ we have three degenerate pions which 
reproduce the masses computed in when the vacuum is trivial. More precisely we 
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reproduce the formula (4.18) of ^H]; when C2 < (i.e. no Aoki phase) and the formula 
(4.14b) of [13], when C2 > and |e| > 1 (i.e. Aoki phase is possible, but not yet reached). 
When C2 > and |e| < 1, our formulas (jl9l2(J|) do not make sense, because they predict 
negative masses. 

However we can do better. We can approach the critical point trying to follow the vacuum 
(fT^ along a path that merges into the Aoki vacuum. This is obtained when r/ — > 1 faster 
than uj —>■ TT. Therefore, instead of (jl5|) . we use the parameterization: 

1] = 1 + 'X^5rj, a; = TT + x^^- (21) 

In this way we get for the pion masses: 

^ = X'4^\ ^ = X^(V^- 16(L86 - W,, + Z,,)). (22) 

Now, setting ^cu 0, we find the expected massles charged pions. Instead, the neutral pion 
is not as expected. The value of m^o in ()22j) agrees with (4.14a) in |13j only for vanishing 
e. Also missing is a better understanding of the transition described in figure [71 over which 
we have no analytical control. On the whole, however, the formula ()17|) provides more 
informations about the critical region, than one could expect. 



5 Conclusions 

We have proposed a method to organize a ChPT calculation which is convenient in the case 
of tmQCD and has a better range of applicability than the ordinary expansion around the 
trivial vacuum. We have computed the vacuum orientation to NLO and compared with 
numerical results which give a rough estimate of the applicability of the expansion. We have 
checked that - within the range of validity of a ChPT expansion - no other orientations are 
possible for the vacuum, apart from combinations of I2 and t^. As long as we are not too 
close to the point 77 = 1 and a; = tt (i.e. where Bqui and W^a are equal in size and point in the 
same direction), the potential has a single minimum, which is essentially the LO vacuum plus 
corrections. When we approach the critical region, a more complicated structure of minima 
is possible, and any perturbative description breaks down. We show - mainly numerically - 
as the LO vacuum is smoothly connected to an Aoki phase, if present. If no Aoki phase is 
possible, we displayed two different scenarios (one smooth, and one involving a transition) 
in which the LO vacuum at non zero twist is connected to the trivial vacuum. 

In the second part we used the vacuum computed before to determine the pion masses. 
We noticed that the introduction of a new degree of freedom allows in principle a better 
determination of the LEC's. Having performed the computation to O(a^), we could also 
determine the splitting between the pion masses in the twisted (m^o) and in the untwisted 
(m^±) directions. Surprisingly, this splitting - which should be measurable even for quite 
heavy pion masses - gives direct information about the phase structure of the system in the 
deep unphysical region dominated by lattice artifacts. Finally we tryed to extrapolate our 
formulas for the pion masses inside the critical region, obtaining some partial results. 
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